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ftjgTHOD and Apparatus for 
mic Sign al Processing and Exploration 



Jf chnical Field 

5 This invention relates to the general subject of seismic exploration and, in 

particular, to apparatus and methods for. the exploration and production of oil and 
gas by identifying structural and stratigraphic features in three dimensions. 



RaHfgrotinrt Of The Invention 

10 In seismic exploration, seismic data is acquired along lines that consist of 

geophone arrays onshore or hydrophone streamer traverses offshore. Geophones 
and hydrophones act as sensors to receive energy that is transmitted into the 
ground and reflected back to the surface from subsurface rock interfaces. Energy 
is usually provided onshore by Vibroseis® vehicles which transmit pulses by 

15 shaking the ground at pre-determined intervals and frequencies on the surface. 
Offshore, airgun sources are usually used. Subtle changes in the energy returned 
to surface often reflect variations in the stratigraphic. structural and fluid contents of 
the reservoirs. 

In performing three-dimensional (3-D) seismic exploration, the principle is 
20 similar, however, lines and arrays are more closely spaced to provide more detailed 
sub-surface coverage. With this high density coverage, extremely large volumes of 
digital data need to be recorded, stored and processed before final interpretation 
can be made. Processing requires extensive computer resources and complex 
software to enhance the signal received from the subsurface and to mute 
25 accompanying noise which masks the signal. 

Once the data is processed, geophysical staff compile and interpret the 3-D 
seismic information in the form of a 3-D data cube (See FIG. 1) which effectively 
represents a display of subsurface features. Using this data cube, information can 
be displayed in various forms. Horizontal time slice maps can be made at selected 
30 depths (See FIG 2). Using a computer workstation an interpreter can also slice 
through the field to investigate reservoir issues at different seismic horizons. 
Vertical slices or sections can also be made in any direction using seismic or welt 
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data. Time maps can be converted to depth to provide a structural interpretation at 
a specific level. 

Seismic data has been traditionally acquired and processed for the purpose 
of imaging seismic reflections. However, changes in stratigraphy are often difficult 
5 to detect on traditional seismic displays due to the limited amount of information 
that stratigraphic features present in a cross-section view. Although such views 
provide an opportunity to see a much larger portion of these features, it is difficult to 
identify fault surfaces within a 3-0 volume where no fault reflections have been 
recorded. 

10 Coherence and sembience (a measure of multichannel coherence) are two 

measures of seismic trace similarity or dissimilarity. As two seismic traces increase 
in coherence, the more they are alike. Assigning a coherence measure on a scale 
from zero to one, a 0" indicates the greatest lack of similarity, while a value of "1" 
indicates total or complete similarity (i.e.. two identical traces). Coherence for more 

1 5 than two traces may be defined in a similar way. 

One method for computing coherence was disclosed in a U.S. Patent 
Application by Bahorich and Farmer (assigned to Amoco Corporation) having a 
serial number of 08/353,934 and a filing date of December 12. 1994. A method for 
computing semblance was disclosed in a U.S. Patent Application by Marfurt et al. 
20 (assigned to Amoco Corporation) having an application number of 60/005.032 and 
a filing date of October 6, 1995. The Marfurt et at. invention included a brute force 
search over candidate dips and azimuths. 

As good as both methods have proved to be, they have some limitations. 
Improved resolution and computational speed are always desirable. 

25 

Summary O f The Invention 

In accordance with the present invention, a multitrace eigen-decomposition 
process is disclosed that is more robust and that has higher resolution than 
previously known methods. In one embodiment of the invention, a method is 
30 disclosed for the exploration of gas and oil. The method comprises the steps of 
accessing a data set of seismic signal traces distributed over a pre-determined 
three-dimensional volume of the earth; in a running window determining the outer- 
product of at least two data vectors formed from at least two seismic traces; 
forming a covariance matrix by adding the outer-products; computing a seismic 
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attribute that is a functiorr of at least the dominant eigenvalue of the covariance 
matrix; and forming a map from the computed seismic attributes over at least a part 
of the pre-det rmined three dimensional volum of the earth. 

In another embodiment of the invention, the process of the invention is 
5 encoded on a computer readable media (e.g., magnetic disk, magnetic tape, 
CD-ROM, etc.) for directing the operation of a computer to compute the seismic 
attributes. In other embodiments of the invention, a map is prepared from the above- 
described process and the map is used to locate oil and gas deposits. 

This technique is particularly well suited for interpreting fault planes within a 
10 3-0 seismic volume and for detecting subtle strattgraphic features in 3-D. This is 
because seismic traces cut by a fault line generally have a different seismic 
character than seismic traces on either side of the fault Measuring seismic trace 
similarity, (i.e., coherence or 3-D continuity) along a time slice reveals lineaments of 
low coherence along these fault lines. Such coherency values can reveal critical 
15 subsurface details that are not readily apparent on traditional seismic sections. 
Also by calculating coherence along a series of time slices, these fault lineaments 
identify fault planes or surfaces. 

Numerous other advantages and features of the present invention will 
become readily apparent from the following detailed description of the invention, the 
20 embodiments described therein, from the claims, and from the accompanying 
drawings. 



Brief Description of the Drawings 

FIG. 1 is a pictorial representation of the information obtained from 
25 processing 3-D seismic data; 

FIG. 2 is a pictorial representation of a horizontal time slice of 3-D seismic 
data processed in accordance with the prior art 

FIG. 3 is a drawing depicting two adjacent seismic traces; 

FIG's. 4 through 8 are schematic diagrams depicting the coherency of a 
30 pair of seismic traces in accordance with the present invention; 

FIG. 9 is a pictorial representation of a running window analysis cube; 

FIGs. 10A, 10B and 10C are schematic diagrams of groups of seismic 

traces; 
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FIG's. 11 A, 11B and 11 C are schematic diagrams of two dimensional 
analysis windows; 

FIG. 12 is an elementary process flow diagram; and 

FIG's. 13A, 13B and 13C are pictorial representations of the same 
5 horizontal time slice in accordance with the inventions of Bahorich et al., Marfurt et 
al. and the present invention. 



Detailed Description 

White this invention is susceptible of embodiment in many different forms, 
10 there is shown in the drawings, and will herein be described in detail, several 
specific embodiments of the invention. It should be understood, however, that the 
present disclosure is to be considered an exemplification of the principles of the 
invention and is not intended to limit the invention to the specific embodiments or 
algorithm so described. 

15 Before describing the process of the invention, the underlying principle of 

operation will be described. 

Consider two traces t, and tj over a specified time window or depth 
window of N samples for which coherence is to be evaluated. A representative 
diagram of the traces and the relevant analysis window is shown in FIG. 3. The 

20 first trace t, consists of the time series (t 11t t 12 t IM ) and the second trace tt 

consists of the time series (t 21 , t^). In these two time series, the first index 

refers to the trace number (i.e., trace 1 or trace 2), while the second index refers to 
the sample number. 

By plotting one trace against the other in the familiar two-dimensional 
25 Cartesian coordinate system a better understanding of the meaning of coherence in 
the context of the present invention can be obtained. Plotting equivalent time 
samples from the two traces [i.e.. the point pairs (t„, t 21 ), (t 12 , t^), .... (t 1Nl t^)], 
produces a crossptot of the two time series. Letting the x-axis represent the first 
trace t, and the y-axis the second trace results in the diagram depicted in FIG. 
30 4. It is the pattern formed by these points which manifests the coherence of the two 
traces. 

The general shap f these two correlated traces is a set of points 
represented by an Hips . This ellipse is a generalization since it does not 



BNSDOCID: <WO 9739367A1_L> 



WO 97/39367 - PCT/US97/00249 

represent each individual point but rather describes the "overall" nature of all the 
points. The major and minor axes of this ellipse will be oriented in a direction which 
is determined by the g metry of the paired points. The lengths of the two axes are 
also determined by this geometry. A typical representation of these points and the 
5 corresponding ellipse is shown in FIG. 5. 

The directions and magnitudes of the major and minor axes of the ellipse 
may be represented by two scaled vectors with the longer vector oriented along the 
major axis and the shorter vector along the minor axis. The magnitudes of these 
two vectors correspond to the two eigenvalues of the data covariance matrix and 
10 the normalized vectors correspond to the eigenvectors. The eigenvectors, scaled 
by their respective eigenvalues, denote the magnitudes and directions of the major 
and minor axes. The •principal component" corresponds to eigenvector which is 
associated with the dominant eigenvalue. 

The figures that follow (FIGs. 6 through 8) are intended to give an intuitive 
15 understanding of the mechanics behind the previous discussion. In these figures, 
traces are constructed using simple amplitude and phase variations, and the effect 
of these variations is observed on the associated eigenvalues and eigenvectors. 

FIG. 6 demonstrates how two identical traces degenerate to a 45 degree 
line. (i.e. t an ellipse with the minor axis collapsed to zero). The "zero length" 
20 indicates that the second eigenvalue is "zero" and indicates that the eigenvector 
corresponding to the dominant eigenvalue is aligned with the major axis. The 
coherence is maximum, with a value of one. 

The situation for two traces having equal amplitudes and a 45 degree 
phase difference is shown in FIG. 7. This demonstrates how a phase shift 
25 lengthens the minor axis and therefore increases the magnitude of the second 
eigenvalue. The two eigenvectors scaled by their respective eigenvalues are also 
shown. Due to the difference in these traces, the coherence is reduced to a value 
of less than one. 

Finally, in FIG. 8 both the phase and amplitude are allowed to vary. The 
30 two traces have both a 45 degree phase shift and a 2 to1 amplitude ratio. The 
resulting ellipse has a nonzero minor axis (the second eigenvalue is nonzero), 
reflecting the phase difference. In addition, the ellipse and eigenvectors are rotated 
due to the amplitude difference. Again, the amplitude and phase variations produce 
a reduction in coherence. 
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The main point of th forgoing discussion is to show heuristicalty that 
coherence can be expressed as a function of the eigenvalues, X 1 and \ 2 , and the 
eigenvectors v, and v 2 . Functi nally, an expression for coherence is: 

Coherence =* / (X,. X 2 , v 1t v 2 ). (1) 

5 This procedure for two traces can easily be extended to any desired 

number of traces. From a practical standpoint the computational burden increases 
with an increasing number of traces and limitations are imposed onty by the 
available computing power. For a 3-D seismic data set this analysis can be 
repeated over a moving or running analysis window in space and time (or space 
10 and depth), resulting in a measure of coherence at the center of the moving 
window. The output is a 3-D data set consisting of coherence values defined over 
the original data volume. 

It will become apparent from the following discussion that one advantage of 
this process is that different aspects of the data are distributed among the 
15 eigenvalues and the eigenvectors. Information, such as amplitude and phase, can 
now be analyzed and dealt with in a robust and rigorous way. Resolution, for 
example, can be improved by manipulating the eigenvalues and eigenvectors. The 
benefits can be observed visually in the computed coherence slices. Another 
important aspect of coherence values in accordance with the present invention is 
20 that they exhibit sensitivity to discontinuities in the original data and reveal subtle 
geological features, such as faults and channels. 

Returning to the process of the present invention, the first step is to obtain 
or access a set of three dimensional seismic data. Such data is in the form of 
seismic signal traces distributed over a three dimensional volume of the earth. 
25 Methods by which such data is obtained and reduced to digital form for processing 
as 3-D seismic data are well known to those skilled in the art Such data is 
routinely acquired by geophysical vendors who specialize in land surveys or ocean 
surveys. Such data is also sold or licensed by vendors and is generally conveyed 
or stored on magnetic tape for transfer into the memory of a seismic work station. 

30 The next step is to divide the 3-D data set into a plurality of cells or analysis 

cubes 20 (See FIG. 9). These cubes 20 perform the function of dividing or sorting 
the seismic data into groups or cells for further processing. In effect one analysts 
cube sweeps through th entire seismic data set or 3D data cube 30. Each 
analysis cube 20 comprises a stack of generally flat rectangular time layers 22. 
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For the purpose of simplicity, the 3-D data cube 30, the analysis cube 20. 
and the tim layers 22 are shown in the form of right-angled parallelepipeds or 
cubes (g n rically a "ceil'). Those skilled in the art will appreciat the simplicity of 
rectangular geometry in performing repetitive operations over 3-D data set. Other 
5 geometry and cell shapes are possible and may be warranted under the 
circumstances. 

Returning to FIG. 10, each time layer 22 or slab has portions of seismic 
traces t, (only one being shown to avoid cluttering the drawings) passing 
therethrough. Nine traces, laid out in a 3 by 3 evenly spaced grid are depicted in 
10 FIG. 9. Five traces in a star pattern (See FIG. 10A) or three traces (See FIG. 10B) 
may be used. A symmetric arrangement is preferred. Nine traces in each time 
layer are often better than three traces. Two traces are the least that can be used. 

To help visualize the concept of the invention, the reader is refered to 
FIG.s 11A and 11B. In particular, it is often difficult to visualize these concepts in 
15 nine dimensions (as is the case of 9 traces). In FIG.s 11A and 11B. a two 
dimensional running window 24 (or analysis window) is shown with just two traces 
t, and ta contained therein. Each trace (See FIG. <1C) comprises a time series of 
N samples 

20 t-(t t ,.t 99i ....W N ) 

The analysis window 24 of FIG. 11C is further divided into a plurality of vertically 
stacked, rectangular time layers 22. 

Before proceeding, it should be understood that in choosing window size 
and spacing, there is a tradeoff between resolution and stability, in other words. 
25 smafl analysis windows or cubes allow higher spatial or temporal frequency in the 
resulting parameter estimate, but give less statistical stability or fewer degrees of 
freedom to those estimates. On the other hand, very large windows have poor 
resolution and tend to smear the data such that important geologic features may be 
lost 

30 Returning to FIG. 11C. the data points within each time layer 22 define 

seismic data vectors (here a 1 by 2 matrix, where N=2). As such, the two 
dimensional time layers of FIG. 11C (or the three dimensional time layer of FIG. 9) 
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form or define vector windows. Thus, nine traces (N=9) will result in a data vector 
having nine elem nts. 

The outer product of the data vector within each vector window or time 
layer 22 results in a N by N matrix. So a nine element data vector gives rise to a 9 
5 by 9 matrix. Adding these matricies (one matrix for each vector window) results in 
a N by N covariance matrix for the entire analysis window 24 (or the analysis cube 
20). So if the analysis cube comprises nine time layers, nine N by N matricies are 
added together to form one N by N covariance matrix. 

In accordance with the present invention, a very useful and non-obvious 
10 measure of trace coherency is obtained by computing the eigenvalues of the 
covariance matrix. In particular, the largest or dominant eigenvalue of the 
covarience matrix and the sum of the eigenvalues of the covariance matrix are 
computed. The ratio of these two numbers represents the size of the dominant 
eigenvalue relative to the sum of eigenvalues. It also indicates the variability of 
15 traces within the analysis cube. Expressed mathematically, a useful seismic 
attribute is represented by: 

^-=4* (2) 

t - I 

where X, are the eigenvalues of the covariance matrix, and X, is the dominant 
eigenvalue. As such, ¥ is an indication of how well the point elements of the 
20 seismic data vectors are correlated (See FIGs. 4 and 5). 

Each and every eigenvalue of the respective covariance matrix need not be 
expressly calculated. Those skilled in the art know that methods exist for 
calculating just the dominant eigenvalue (e.g., the power method, the Rayleigh 
quotient [the faster of the two], etc.). Moreover, it is also known that the sum of the 
25 diagonal elements of the covariance matrix is equal to the sum of the eigenvalues 
of that covariance matrix. For convenience, the ratio (i.e., a measure of coherency) 
of equation (2) can be assigned to the center of the analysis cube 20 or the 
analysis window 22. It should be appreciated that one advantage of using the 
dominant eigenvalue is that it tends to more directly show the variability of the 
30 traces within the analysis window. The dominant eigenvector would not be a 
measure of the variability (i.e., the coherence). 
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In performing this measure of coherency, the analysis cube 20 or analysis 
window 22 effectiv ly sweeps (i.e., laterally and vertically) accross the entire 3-D 
volume of FIG. 9 or all the traces of FIG. 1 1 A. Pr ferably, adjac nt analysis cubes 
20 or vector windows 22 overlap one another (See windows 24, 24* and 24* of FIG. 
5 11 A). Overlap improves spatial resolution. 

The result of the analysis cube 20 sweeping over the entire 3D data volume 
30 and the assignment of the coherency measure of equation (2) is an array of 
coherency values assigned along each trace at the location of each data vector. 
As such, the 3D data volume is converted into a 3D "coherency cube." 

10 The data or coherency measures contained within the coherency cube is 

most conveniently interpreted by displaying the coherency data in the form of a 
seismic attribute map. Such a map is most often in the form of a display of those 
coherency values lying along a surface passing through the coherency cube. Two 
examples are a flat plane passing through a common horizontal time slice, and a 

15 curved surface passing through a seismic horizon line selected by a seismic 
interpreter. Another example is a line representative of geologic depositions time 
so as to capture features from the same geologic age. 

The coherency values are easily displayed for interpretation when 
displayed as a shade of gray (e.g., white indicating the highest coherency and black 

20 indicating the lowest coherency) or another color scale. (See US Patent 4,970.699 
for a "Method for Color Mapping Geophysical Data"). 

Landmark and GeoQuest interpretive workstations, for example, can be 
used to view and interpret faults and stratigraphic features by loading the 
coherency cube as a seismic volume. Such workstations are commonly used by 

25 those skilled in the art Unprocessed 3D seismic data can be conveniently loaded 
onto the workstation by means of a magnetic tape or disk which is encoded with 
instructions for the computer to perform the above-described process. 
Visualization software (e.g., Landmarks's SeisCube software) may be employed to 
rapidly slice through the coherency cube to aid in understanding complex fault 

30 relationships. Coherency displays, including printouts in the form of seismic 
attribute maps, can reduce interpretation cycle time when used in-selecting which 
seismic lines to interpret enabling the interpreter to work around poor data areas. 
In addition, subtle stratigraphic features and complex faulting, which are not readily 
apparent on traditional seismic displays, can be rapidly identified and interpreted. 
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FIG'S. 13A, 13B ^nd 13C provide comparisons of the same seismic 
information displayed and processed by other processes and in accordanc with 
the present invention. The differences are readily apparent. 

Of course, the process of the invention is most conveniently carried out by 
5 writing a computer program to carry out the steps just described. Such processes 
are routinely carried out on work stations previously identified. 

An elementary process flow diagram is illustrated in FIG. 12. In one 
embodiment of the invention, a software program is written in FORTRAN 77 to 
perform the process just described. 3-D seismic data 30 is read into memory. 
1 0 Based on the size and content of the available seismic data, an initialization step 32 
is performed, program parameters are initialized, data ranges are established, 
preliminary checks are made, and window sizes are set. Defaults are read or user 
selected options are read. Next processing 34 is begun. In particular, subroutines 
36 are called to sweep the data volume with an analysis cube. Within each 
15 analysis cube, a subroutine 38 computes the covariance matrix, and another 
subroutine 40 computes the dominant eigenvalues, the sum of the eigenvalues and 
the resultant coherency value. Finally, the results are combined 42 and the 
computed values are stored 44 in the form of a coherency cube. Thereafter, a 
workstation operator can access the coherency cube to display selected portions 
20 (e.g., time slice through the cube) on a CRT 46, to create a printout or seismic 
attribute map 48, to perform further analysis or to transfer to memory or to a tape 
50 for further processing elsewhere. 

Those skilled in the art are cautioned to use the method of the invention 
with due regard to what is known about the stratigraphy and geology of the region 
25 covered by the 3-D survey. 

Coherency maps have been run on several 3-D surveys. At depths of 
reasonable data quality, approximately 90% of the faults can be readily identified. 
Faults were identified on coherency maps which were very subtle on seismic 
sections, but clearly present on the coherency maps because of the robustness of 
30 the method and the map perspective of fault patterns. Since coherency maps can 
be run on uninterpreted time slices, the present invention offers a means to greatly 
accelerate mapping of the structural framework and to reveal details of fault 
relationships which would otherwise be interpreted only through tedious fault 
picking. 
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Coherence maps were generated along picked horizons 
and clearly identified salt/shale diapirs offshore 

In other locations mud and gas volcanos were clearly 
5 indicated using the process of the invention. 

Several coherence time slices showed remarkable details 
of stratigraphic features, such as abandoned river channels, mud 
flows, point bars and submarine canyons. On seismic sections, 
these features were sometimes apparent but, in some cases, were 
10 unidentifiable even with close scrutiny. 

This invention, like that of Bahorich et al. and Marfurt et al., provides a 
method of revealing fault planes within a 3-D volume where no fault reflections have 
been recorded. Faults are often critical to the accumulation of oil. A fault may form 
a seal by cutting off a structural or stratigraphic feature so the oil is trapped against 
15 the fault On the other hand, if the fault plane contains rubble that has not been 
cemented, it may form a conduit for fluids. This may allow the hydrocarbons to drift 
up the fault plane into the feature and be trapped in it or to escape from the feature 
by drifting up the fault plane out of it Thus, fault tines can predict flow patterns in a 
reservoir and communication between injector and producing wells, for example. 

20 Seismic discontinuities can also provide the needed link to enable reservoir 

prediction between the wells and establish reservoir continuity and flow patterns 
across a field. Coherency technology can be used for finding, identifying and 
mapping of subsurface structural and sedimentological features such as faults, salt 
diapirs, unconformities, channel systems, karsting and carbonate reef fades which 

25 are commonly associated with the entrapment and storage of hydrocarbons. 
Therefore, this technology aids in the finding, extraction and production of 
hydrocarbons. In addition, it is used for identifying both shallow and deep drilling 
hazards (e.g., places where there is gas that is too close to the surface or where 
there are instabilities). Still another example is the use of the invention to search 

30 for leakage paths from known reservoirs or underground storage caverns. 

Coherency mapping of 3-D seismic is an extremely powerful and efficient 
tool for mapping both structure and stratigraphy. The new method is particularly 
sensitive to any lateral variation in wavelet character and therefore is particularly 
sensitive to the common causes of lateral variations in th wavelet (i.e., fault 
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displacement or stratigraphic variations). Thus, the subject invention encompasses 
a process, devices on which the process is recorded in the form of computer 
instructions, the product (e.g.. a map) of that process, and the manner in which 
such a product is used in the exploration of gas and oil. 

5 From the foregoing description, it will be observed that numerous 

variations, alternatives and modifications will be apparent to those skilled in the art. 
Accordingly, this description is to be construed as illustrative only and is for the 
purpose of teaching those skilled in the art the manner of carrying out the invention. 
For example, the seismic traces have been depicted as having equal spacing. 
10 Unevenly spaced traces (See FIG. 10C) may be conveniently converted into even 
spacing by interpretation. As another example, it may be useful to filter the traces 
forming the data vector to eliminate outlyers. A median filter may be used to sort 
the elements of each data vector. Edge cutoffs may be used to achieve further 
smoothing. 

15 Moreover, other algorithms may be used to measure the similarity of 

nearby regions of seismic data or to generate the "coherence cube." The 
coherency value or seismic attribute of equation (2) serves as a rather robust 
estimate or measure of signal discontinuity within geologic formations as well as 
signal discontinuities across faults and erosional unconformities. Other 

20 combinations of the eigenvalues of the covariance matrix are suggested (e.g., 
arithmetic mean, root mean square, average, median, square root of the sum of the 
squares, square root of the product of the squares, minimum, maximum, sum, 
product etc.). Moreover, the process of the invention may be combined with other 
attributes (e.g.. AVO slope, etc.) and also applied to multi-component seismic data. 

25 Also certain features of the invention may be used independently of other 

features of the invention. For example, geologic features identified in accordance 
with the present invention can be overiayed with a velocity map to provide a means 
of cross checking the velocities. 

In addition, white coherence slice maps by themselves are very powerful 

30 mapping toots, when used in conjunction with reconnaissance mapping of 
amplitudes and dip maps, there is promise of a technological milestone in mapping 
effectiveness for the Gulf of Mexico or similar basins with 3-D seismic data readily 
available. It is believed that detailed mapping of structure and stratigraphy yvill be 
accelerated by mapping in a map view and less by traditional line by tine picking. 
35 Interpretati n in a map view f "reconnaissance" data offers significant 
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improvement in quality and quantity of interpretation. Moreover, the process of the 
invention is inherently fast Such speed aids in rapidly making bidding choices 
when concessions becom available. 

Finally, it should also be understood that the principle of the invention 
5 should be equally applicable to other fields (e.g., passive sonar, in which case the 
sensors would be acoustic and the signal sources could be hostile submarines; 
earthquake and nuclear weapon detonation detection systems, in which case the 
sensors would be seismic and the signal sources could be earthquake or explosion 
epicenters; astronomical interferometry, wherein the sensors would be 
10 radiotelescopes and the signal sources could be distant galaxies or quasars; and 
phased-array radars in which case the sensors would be the array antennae) 
wherein signals (e.g., radar, sonar, radio frequency energy, etc.) are processed to 
form images or locate changes in the structure represented by such images. 

Thus, it will be appreciated that various modifications, alternatives. 
15 variations, and changes may be made without departing from the spirit and scope 
of the invention as defined in the appended claims. It is, of course, intended to 
cover by the appended claims all such modifications involved within the scope of 
the claims. 
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CLAIMS 

We claim: 

1 . A method for the exploration of hydrocarbons, comprising the steps 

of: 

5 a) obtaining a set of seismic traces distributed over a pre- 

determined three-dimensional volume of the earth; 

b) dividing said three-dimensional volume into a plurality of 
vertically stacked and generally spaced apart horizontal time layers and arranging 
said time layers into a plurality of cells that extend laterally and vertically, each of 

10 said time layers having portions of at least two seismic traces located therein that 
define a data vector 

c) calculating in each of said time layers of said cells the 
outer-product of said data vector 

d) combining said outer products to obtain a covariance 
1 5 matrix for each of said ceils; 

e) computing in each of said ceils a measure of the 
coherency of said seismic traces, wherein said measure of coherency is at least a 
function of the largest eigenvalue of said covariance matrix; and 

f) forming a seismic attribute map from a plurality of said 
20 measures of coherency of said seismic traces. 

2. The method of Claim 1, where in step (f) said map is formed by 
displaying said measures of coherency relative to a surface passing through a pre- 
determined seismic horizon. 

25 

3. The method of Claim 1, where in step (f) said map is formed by 
displaying said measures of coherency relative to a surface passing through a pre- 
determined time line. 

30 4. The method of Claim 1. where in step (b) said cells comprise 

analysis cubes having at least five seismic traces located th rein; and wh re in 
performing step (c) each outer product is in the form of a 5 by 5 matrix. 
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5. The method of Claim 4, where in step (b) said cells comprise 
analysis cubes that hav portions of at least nin seismic traces located th rein; 
and wherein said data vectors have nine elements. 

5 

6. The method of Claim 5, where in step (b) said nine seismic traces 
are arranged into a three by three grid. 

7. The method of Claim 1, where in step (b) said cells are less than 
10 100 milliseconds in thickness. 

8. The method of Claim 1 , wherein step (c) is performed in the time 

domain. 

15 9. The method of Claim 1. where in step (e) is performed by: 

computing said largest eigenvalue of said covariance matrix, computing the sum of 
the eigenvalues of said covariance matrix, and computing the ratio of said largest 
eigenvalue to said sum of the eigenvalues of said covariance matrix. 

20 10. The method of Claim 9, wherein said sum of the eigenvalues of 

said covariance matrix is computed by forming the sum of the diagonal elements of 
said covariance matrix. 

11. The method of Claim 1, where in performing step (b) one of said 
25 two seismic traces in each cell is located in an adjacent cell, such that said cells 
spatially overlap one another. 



WO 97/39367 

16 

12. A method of locating subterranean features, faults, and contours, 
comprising th steps of 

a) acquiring 3-D seismic data covering a pre-determined 
volume of the earth, said data comprising seismic traces characterized by time, 

5 position and amplitude; 

b) dividing at least a portion of said volume into at least one 
array of relatively small, adjacent overlapping, three-dimensional analysis cubes, 
wherein each of said analysis cubes contains at least three laterally separated 
seismic traces, and dividing each analysis cube into a plurality of sample intervals 

10 such that each sample interval defines a plurality of one by three data vectors; 

c) computing a seismic attribute for each said cube that is a 
function of the dominant eigenvalue of a covariance matrix formed from outer 
products of said data vectors; and 

d) storing said seismic attributes of said analysis cubes for 
1 5 display in the form of a two-dimensional map of subterranean features. 

13. The method of Claim 12, where in step (c) said seismic attribute is 
a function of the ratio of said dominant eigenvalue to the sum of at least two of the 
eigenvalues of said covariance matrix of said cube. 

20 

14. The method of Claim 13, where in step (c) said seismic attribute is 
a function of the ratio of said dominant eigenvalue to the sum of all the diagonal 
elements of said covariance matrix. 

25 15. The method of Claim 14, where said seismic attribute is assigned 

to the center of its analysis cube. 

16. The method of Claim 15. wherein step (b) is performed on a 
plurality of time slices further including the step of: 
30 e) displaying said seismic attributes of successive time slices 

passing through said centers of said analysis cubes, to identify relativ space and 
tim invariant features. 
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17. In seismic exploration wherein 3-D seismic data comprising 
reflected s ismic en rgy is recorded as a function of time to produce a sen s of 
seismic traces, and wherein a computer is used that is adapted to process such 
5 seismic traces, an article of manufacture comprising: 

a medium that is readable by a computer and that carries instructions for 
said computer to perform a process comprising the steps of. 

(a) accessing 3-D seismic data over a predetermined volume, 
said data comprising seismic signal vectors characterized by time, position and 

10 amplitude; and 

(b) ascertaining the similarity of nearby regions of said 3-D 
seismic data of said volume by: 

(1) dividing at least a portion of said data into an array 
of relatively small, adjacent, overlapping, three-dimensional analysis cubes. 

15 wherein each of said analysis cubes contains at least two data vectors; and 

(2) computing a seismic attribute for each cube that is 
a function of the principal eigenvalue of a covariance matrix that is formed from a 
sum of outer products of said vectors of said cube. 

20 18. The article of manufacture of claim 17. wherein said medium 

carries instructions for said computer to perform step (2) by computing the ratio of 
said principal eigenvalue to the sum of eigenvalues of said covariance matrix. 

19. The article of manufacture of claim 17. wherein said medium 
25 carries instructions for said computer to perform step (2) by computing the ratio of 

said principal eigenvalue to the sum of the diagonal elements of said covariance 
matrix. 

20. The article of manufacture of claim 19, wherein said medium 
30 carries instructions for said computer to perform step (1 ) by forming analysis cubes 

having a generally rectangular array of least five seismic traces located therein; and 
wherein said covariance matrix is at least a five by fiv matrix and is formed from at 
least three outer product matrices. 
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21. The article of manufacture of claim 20, wherein said medium 
carries instructions for said computer to assign said seismic attributes to the center 
of its analysis cube. 

5 

22. In seismic exploration wherein reflected seismic energy is recorded 
as a function of time to produce a series of seismic traces, a method comprising the 
steps of 

(a) determining the outer-product two data vectors formed 
1 0 from at least two seismic traces; 

(b) forming a covariance matrix by adding said outer-products 

of step (a); 

(c) computing a seismic attribute that is a function of at least 
the principal eigenvalue of said covariance matrix of step (b); 

15 (d) repeating steps (a) through (c) across at least a part of at 

least one time window; and 

(e) forming a map of said seismic attributes over said time 

window. 



20 23. The method of Claim 22, wherein step (c) is performed by 

computing the ratio of said principal eigenvalue to at least a partial sum of the 
eigenvalues of said covariance matrix. 

24. The method of Claim 22, wherein step (c) is performed by 
25 computing the ratio of said principal eigenvalue to at least a partial sum of the 
diagonal elements of said covariance matrix. 



25. The method of Claim 22. wherein step (d) is performed by using at 
least one seismic trace from the previous performance of step (a) and at least two 
30 new seismic traces that are located adjacent to said at least one seismic trace. 



BNSDOCID: <WO 9739367A1J_> 



WO 97/39367 




PCT/US97/00249 



19 

26. The method of Claim 22, wherein st p (a) comprises the steps of 

(1) accessing 3-D seismic data over a predetermined volume 
of the earth, said 3-D seismic data s ismic traces characterized by time, position 
and amplitude; and 

5 (2) dividing a portion of said volume into at least one time 

window comprising an array of relatively small, overlapping, three-dimensional 
analysis cubes that contain said at least two seismic traces. 

27. A method of seismic exploration, comprising the steps of: 

10 a) reading a 3-D seismic data set comprising seismic signal 

traces distributed over a volume of the earth; 

b) selecting at least one horizon slice from said volume and 
forming thereon cells that are arranged into laterally extending rows and columns, 
each of said cells comprising at least three seismic traces extending generally 

15 therethrough; 

c) computing for each of said cells: 

(1) the outer-product of data vectors defined by a 
plurality of time intervals on each side of the center said cell; 

(2) a covariance matrix from said outer products of 

20 step (1); and 

(3) at least the largest eigenvalue of said covariance 

matrix; and 

(d) examining said eigenvalues of said cells along said at least 
one horizon slice. 

25 

28. The method of Claim 27, wherein step (3) is performed by 
displaying a representation of said largest eigenvalues of said ceils across at least 
one horizontal time slice. 

30 29. The method of Claim 28, wherein said representation is a function 

of said largest etgenvalu of said cell and the sum of the eigenvalues of said 
covariance matrix of said cell. 
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30. In seismic exploration wherein reflected seismic energy is recorded 
as a function of time to produce a series of seismic trac s. a method comprising the 
steps of. 

5 (a) arranging 3-D seismic data into relatively small, 

overlapping, three-dimensional analysis cubes that contain a plurality of seismic 
traces; 

(b) determining the outer-product of the data vectors defined 
by said analysis cubes; 

10 (c) forming a covariance matrix for each cube by adding said 

outer-products of step {b); 

(d) computing a seismic attribute that is a function of the ratio 
of the principal eigenvalue of each covariance matrix to the sum of all the 
eigenvalues of that covariance matrix; and 
15 (e) arranging said seismic attributes for being displayed as a 

map. 

31 . A device, comprising: 

computer readable means carrying instructions for a process comprising the 
20 steps of: 

(1) reading into memory 3-D seismic data that covers a pre* 
determined volume of the earth; 

(2) digitally sorting said 3-D seismic data into an array of 
relatively small three-dimensional cells wherein each of said cells contains at least 

25 three seismic traces; 

(3) calculating in each said cells a coherency value from the 
eigenvalues of a covariance matrix formed from a plurality of outer products of said at 
least three traces; and 

(4) storing said coherency values of said cells for displaying 
30 same in the form of two-dimensional map of subterranean features represented by 

said coherency values. 
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32. The device of Claim 31, where in step (3) said coherency value is at 
least a functi n of the largest of said eigenvalues of said covariance matrix. 

33. The device of Claim 32, wherein said coherency value is a function of 
5 said largest eigenvalue and a sum of said eigenvalues. 

34. The device of Claim 31 wherein said computer-readable means is 
selected from the group consisting of a magnetic tape, a magnetic disk, an optical disk 
and a CD-ROM. 

10 

35. A method of locating subterranean features, faults, and contours, 
comprising the steps of: 

a) obtaining seismic data covering a predetermined volume of 

the earth; 

15 b) dividing said volume into an array of relatively small three- 

dimensional cells wherein each of said ceils is characterized by at least two seismic 
data vectors located therein; 

c) computing a covariance matrix from the outer products of 
said data vectors 

20 d) mapping a representation of the eigenvalues of said 

covariance matrix. 

36. The method of Claim 35 where in step (c) is performed by using a 
covariance matrix formed from adding a plurality of outer products. 

25 

37. The method of Claim 35, wherein step (d) is performed by mapping 
the ratio of the largest eigenvalue to a sum of eigenvalues. 
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38. A method of prospecting for hydrocarbon deposits, comprising the 
steps of. 

a) obtaining a seismic attribute map of coherency values of 3-D seismic 
data for a predetermined three-dimensional volume of the earth, said map being 
5 generated by using a computer and a program for said computer that instructs said 
computer to perform the following steps: 

(1) reading said data and sorting said volume into an 
array of relatively small three-dimensional cells, wherein each of said cells has at least 
two seismic data vectors located therein; and 

10 (2) calculating in each of said ceils a coherency value for 

said seismic traces that is a function of the eigenvalues of a covariance matrix formed 
from outer products of said data vectors; and 

(b) using said map to identify subsurface structural and sedimentological 
features commonly associated with the entrapment and storage of hydrocarbons. 

15 

39. The method of Claim 38. further including the step of using said map 
to identify drilling hazards. 



40. The method of Claim 39, further including the step of drilling at a 
20 location identified in step (b). 

41. The method of Claim 38, wherein said program instructs said 
computer to perform step (a)(2) by: 

(i) computing the largest eigenvalue of each covariance matrix 
25 and the sum of the eigenvalues of said covariance matrix; and 

(it) computing the ratio of said largest eigenvalue to said sum. 

42. The method of Claim 41, where in performing step (i) said program 
instructs said computer to compute said sum of eigenvalues by calculating the sum of 

30 the diagonal elements of said covariance matrix. 
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43. A seismic^map prepared by a process, comprising the steps of 

(1) accessing, by means of a computer, a dataset comprising 
seismic signal traces distributed over a pre-d t rmined three-dimensional volum of 
the earth; 



that are arranged in time and space, each of said cells having located therein a 
plurality of data vectors; 

(3) computing in each cell a plurality of outer products formed 
from data vectors located therein; 

(4) combining said outer products to form a matrix for each cell; 

(5) computing the dominant eigenvalue of said matrix and the 
sum of the diagonal elements of said matrix; and 

(6) displaying said largest eigenvalue relative to said sum for 
each matrix of a pre-determined group of cells passing through a pre-determined 
surface. 

44. The seismic map of Claim 43, wherein step (6) is performed by 
obtaining the ratio of said dominant value to said sum for each said matrix of said pre- 
determined group of ceils passing through said pre-determined surface. 

45. The seismic map of Claim 43 where in step (2) each of said data 
vectors has at least three elements. 

46. The seismic map of Claim 43, wherein step (4) is performed by 
adding together said outer products. 

47. A map for the exploration of oil and gas, comprising: 

a) a generally flat medium for recording visually perceptable images 
thereon; and 



(2) 



dividing said three-dimensional volume into a plurality of cells 



WO 97/39367 



9 



PCT/US97/00249 



b) a plurality- of images on said medium that are a function of the 
dominant eigenvalue of a covariance matrix that is formed from the outer products of 
a running window of data vectors representing a 3-D seismic survey. 

5 48. The map of Claim 47, wherein said images are a function of the ratio 

of said dominant eigenvalue to the sum of the eigenvalues of said covariance matrix. 

49. The map of Claim 47, wherein said medium is the face of a cathaode 
ray tube. 

50. The map of Claim 47, wherein said images are a function of the ratio 
of said document eigenvalue to the sum of the diagonal elements of said covariance 
matrix. 

15 si. The map of Claim 47, wherein each running window comprises an 

analysis cube comprising a stack of at least three time layers; wherein each time layer 
contains a data vector therein; and wherein said data vector comprises at least three 
seismic trace elements. 

20 52. The map of Claim 51 , wherein said eigenvalues are assigned to the 

center of each analysis cube. 

53. An exploration map formed by a process comprising the steps of: 

a) forming a coherency cube from data vectors of 3-D seismic data, said 
25 coherency cube comprising a three dimensional array of coherency values that are at 

least a function of the dominant eigenvalues of covariance matricies of said data 
vectors; and 

b) displaying said coherency values as an image on a surface in 
accordance with a predetermined transfer criteria. 

30 
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54. The map of Claim 53, wherein said coherency values are assigned to 
thre dimensional coordinates that are generally coincid nt with elements of said data 
vectors. 

5 55. The map of Claim 53, where in step (b) said surface is a flat plane 

and said predetermined transfer criteria is that said flat plane is generally coincident 
with a time slice through said 3-D seismic data. 

56. The map of Claim 53, where in step (a) each coherency value is at 
10 least a function of said dominant eigenvalue and the sum of the eigenvalues of the 

respective covariance matrix. 

57. A device for use by a computer work station of the type used in the 
exploration of oil and gas, comprising a medium readable by a computer and 

15 containing a representation of a coherency cube, said coherency cube comprising 
measurements of the coherency of 3D seismic data, each of said measurements 
being a function of the eigenvalues of a covariance matrix formed from adding at least 
two outer products of at least two seismic data vectors. 

20 58. The device of Claim 57, wherein said data vectors are characterized 

by space and time coordinates; and wherein said measurements of coherency are 
assigned to said space and time coordinates. 

59. The device of Claim 58, wherein each of said measurements is at 
25 least a function of the dominant eigenvalue of the respective covariance matrix. 

60. The device of Claim 59, wherein each of said measurements is at 
least a function of a sum of eigenvalues. 
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